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I.  Introduction 

Theoretical  studies  are  presented  for  two  different  types  of 
laser- induced  gas-phase  rate  processes.  The  first  (Sec.  II.  A)  is 
collisional  ionization,  where  the  radiation  couples  bound  and 
continuum  electronic  states  of  two  colliding  atoms.  The  theory 
relies  on  classical  mechanics  for  a  description  of  the  nuclear  motion, 
but  treats  the  electronic  degrees  of  freedom  quantum  mechanically. 

The  second  process  (Sec.  II.  B)  is  multiphoton  resonances  in  atom- 
atom  scattering,  where  laser  radiation  induces  transitions  to  quasi¬ 
bound  rovibrational  states  during  the  course  of  the  collision.  The 

collision  process  is  considered  in  the  presence  of  two  simultaneous 
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lasers  (intensity  less  than  1  kW/cm  ) ,  with  all  molecular  degrees 
of  freedom  treated  quantum  mechanically.  The  possibility  of  enhanced 
curve  switching  by  this  process  is  investigated.  Finally,  in  Sec. 

Ill,  the  use  of  laser  radiation  to  excite  electrons  from  the  valence 
band  to  surface  states  of  a  semiconductor  is  considered.  This  can 
lead  to  increased  charge  in  the  surface  region  and  hence  a  strong 
Coulombic  interaction  with  an  adspecies,  which  is  calculated  for  a 
model  system. 

II.  Gas-Phase  Processes 
A.  Collisional  Ionization 

There  are  really  two  problems,  the  structural  and  the  dynamical. 
The  structural  part  is  concerned  with  suitable  choices  of  electronic 
basis  functions,  potential  surfaces,  coupling  matrix  elements  and 
the  like.  These  are  massive  problems  which  we  are  not  equipped  to 
address.  Our  attention  will  be  devoted  to  the  dynamical  part  - 
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the  solution  of  the  equations  of  motion  describing  ionization 

events ,  given  the  raw  structural  ingredients. 

The  simplest  situation  one  can  imagine  -  and  we  will  be 

working  in  this  framework  — — 1  is  an  atom-atom  collision  involving 

the  ejection  of  a  single  electron  in  which  the  different  partial 

wave  events  are  decoupled.  Let  us  just  look  at  the  l-th  partial 

ID 

wave.  The  radial  wave  function,^  ' (R)  (R  =  internuclear  distance), 

IP)  IP) 

can  then  be  regarded  as  a  vector  with  the  components  '(R),i|^  '  (R) )  . 
d  stands  for  the  initial  discrete  channel  and  c  for  a  continuum  of 
channels,  one  for  each  possible  energy  of  the  emitted  electron.  It 
is  apparent  that  0  <  £  <  00 .  Now  the  dynamicist  decides  on  what 

electronic  basis  to  work  with  -  let  us  say  the  one  in  which  motional 

or  nonadiabatic  couplings  between  the  different  channels  are  banned, 
usually  called  the  diabatic  representation.  The  structuralist  then 
hands  him  the  corresponding  set  of  ’matrix  elements’,  consisting  of 
the  discrete  and  ionization  threshold  potential  curves,  W2(R)  and 
W^(R)  respectively,  and  the  bound  continuum  couplings, V£ (R) .  If 
for  the  moment  we  regard  e  as  a  denumerably  infinite  descretized 
set,  the  time- independent  Schrodinger  equation  can  be  cast  in  the 
form 
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an  infinite  set  of  coupled  differential  equations.  We  note  here 
that  Eq.  (1)  is  not  only  good  for  the  description  of  ionization. 

It  can  also  be  looked  at  under  the  guise  of  the  problem  of  colli- 
sionally- induced  light  emission.  Eq.  (1)  is  not  easy  to  solve  in 
general.  Various  ways  have  been  devised  to  get  around  it.  First, 
if  one  is  not  interested  in  a  particular  e  channel  but  only  in  the 
total  decay  rate  from  the  d  channel,  a  complex  potential 

W  —  W_  +  iv  (non-Hermitian  and  hence  non-energy  conserving) 

^  £ 

can  be  introduced  to  describe  the  decay  [1] .  One  then  just  solves 


fa-***- 


+  k2\  *<*>- 


') 


* 


(2) 


Alternatively,  one  can  introduce  a  specific  discretization  scheme  and 
pick  out  only  the  e  of  interest  to  do  the  problem,  or  at  most  a  few 
of  the  nearby  channels,  arguing  that,  on  the  basis  of  the  Franck - 
Condon  principle,  only  nuclear  energy  conserving  ionization  events 
are  significant:  In  this  case  W^+e  is  required  to  be  equal  to,  or 
at  least  approximately  equal  to  W2.  When  one  takes  this  approximation 
seriously,  the  two-state  coupled-channels  discretization  procedure 
[2]  results,  in  which  one  just  solves 


;2  .  2h 


Ve  W1+ej 


-  iUm  +  k2 

R2 


=  0 


(3) 


Variations  of  this  theme  involve  the  inclusion  of  a  few  other  nearby 
channels.  The  above  treatments  are  within  the  framework  of  'exact' 
quantum  mechanics.  They  have  been  amply  worked  over  and  demonstrated 
to  have  been  successful  to  a  certain  extent.  We  will  not  pursue  them 


further  here 
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That  quantum  mechanics  furnishes  the  ultimately  correct  approach 
to  these  problems  is  of  course  not  disputed.  It  is  not  even  claimed 
that  the  so-called  semiclassical  methods  are  necessarily  more  labor 
saving.  Nevertheless,  it  is  interesting  to  se?  what  classical 
mechanics  can  do  for  us,  simply  because,  in  some  sense,  it  is  more 
intuitive.  In  certain  cases  it  is  much  easier  to  think  of  dynamics 
in  terms  of  billiard  balls  moving  on  potential  surfaces  rather  than 
colliding  wave  packets  with  probabilities  of  ending  up  in  different 
channels.  But  the  question  of  exactly  how  far  one  can  push  classical 
mechanics  in  molecular  dynamics  has  never  really  been  settled,  nor 
can  we  pretend  to  be  able  to  supply  an  answer  here. 

The  central  idea  of  semiclassical  methods  involves  the  notion 
of  classical  trajectories  (R  as  a  definite  function  of  t) ,  foregoing 
the  notion  of  deterministic  trajectories  given  a  fixed  set  of  initial 
conditions.  This  effectively  means  that  the  dynamics  in  different 
channels  can  interfere;  but  in  each  it  is  governed  by  a  certain  R(t) . 
The  interference  leads  to  quantum  mechanical  results  otherwise  not 
obtainable  from  straight  classical  mechanics. 

We  will  now  apply  these  ideas  to  the  problem  of  collisional 
ionization  [3] .  First,  instead  of  Eq.  (1) ,  we  can  use  the  time- 
dependent  Schrodinger  equation,  since  R  is  already  regarded  as  a 
function  of  t.  In  so  doing,  the  second-order  coupled  differential 
equations  (in  R)  will  become  first-order  ones  (in  t) .  Then  comes 
a  re-examination  of  the  choice  of  basis.  Do  we  still  want  to  work 
with  the  one  used  in  Eq.  (1)?  We  can,  of  course,  but  we  can  also, 
by  a  transformation  of  basis,  turn  what  is  apparently  a  problem  in 
Eq.  (1)  (that  of  the  infinitely  coupled  nature  of  the  equations)  into 
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an  advantage.  We  can  use  a  basis  set,  characterized  only  by  c  and 
diagonal  in  e,  make  e  continuous,  and  turn  Eq.  (1)  into  a  simple 
differential-integral  equation.  This  is  sometimes  called  a 
Lee-model  approach .  Thus 


d(r) 

zet 


U)  =  c(£) 


=  Y 


becomes,  on  transition  to  the  continuum  limit. 


y^  •  ^de * p  (e 1 ) 


DU>  YU) 

uee'  '  e  '  » 


where  the  D^£,  are  couplings  due  solely  to  the  time-derivative 
operator  ^  =  ft  ^  ,  i.e.,  due  to  motional  or  nonadiabatic  coupling. 
The  diagonal  representation  pointed  to  above  is  then  called  the 
adiabatic  representation,  p(e')  is  a  suitable  density  of  states. 

The  advantages  of  this  approach  are  threefold:  First,  we  can  get 
single-channel  information  for  all  channels  with  one  equation 
[compare  Eq.  (2)1;  second,  no  explicit  function  relating  e  to  R 
need  be  imposed,  as  in  the  Franck-Condon  approximation  e(R)-W2-W1); 
third,  no  truncation  of  the  coupled  discretized  channels  need  be 
used,  in  fact, discretized  channels  are  not  used  at  all. 

(p\  g  g 

Dee'  dePends  explicitly  on  ^(Wj-W^  and  Tt*e  '  an<*  hence  on 
the  dynamics.  Now  two  questions  arise:  What  is  the  potential  that 
one  should  use  for  the  specification  of  R(t);  and  is  it  channel  de¬ 
pendent?  Answers:  in  general,  a  complex  nonlocal  potential  with 

dominant  part  W  (t)  =  W,  (R)  +  +  e;  and  strictly  speaking,  yes. 

1  R* 

The  main  assumption  concerning  R(t)  we  make  is  that,  in  the  solution 

ip\ 

for  a  particular  y '  ’ ,  all  dynamics  is  governed  by  a  single  trajectory 
R  (t)  determined  from  W  (t) . 

£  •  t 
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Eq.  (5)  permits  an  analytic  first-order  perturbative  solution 
that  reveals  clearly  the  contribution  due  to  the  dynamics.  Higher- 
order  solutions  can  in  principle  be  obtained  through  an  iterative 
procedure  but  have  not  been  attempted.  Even  the  evaluation  of  the 
first-order  term  can  be  significantly  simplified  by  the  stationary- 
phase  approximation.  Thus  an  integral  over  t  is  reduced  to  a  sum 
over  stationary  times,  further  delineating  the  part  played  by 
classical  mechanics  in  this  approach.  Finally,  we  point  out  that  if 
certain  approximations  are  made  so  that  the  nonadiabatic  couplings 
D^e;  become  separable  in  e,  the  differential-integral  equation 
(5)  reduces  to  an  integral  equation  —  a  Vol terra  equation  of  the 

second  kind  -  the  exact  solutions  to  which  can  be  much  more  readily 

obtained. 

B.  Multiphoton  Resonances  in  Atom-Atom  Scattering 

In  the  realm  of  radiation-assisted  inelastic  collisions  (RAIC) , 
it  has  become  clear  that  the  exploitation  of  resonances  affords  a 
solution  to  the  problem  of  achieving  good  cross-sections  at  realistic 
laser  intensities [4] .  [By  "realistic"  we  mean  low  enough  to  be 
experimentally  attainable  and  also  below  the  intensity  threshold 
required  for  unwanted  side-effects  like  ionization.]  At  the  same 
time  it  has  also  become  apparent  that  the  non-linear  and  saturation 
regimes,  which  have  traditionally  been  thought  to  require  "intense" 
lasers,  can  become  obvious  at  much  lower  intensities. 

Laser-induced  resonances  involve  transitions  between  pairs  of 
electronic  states,  and  attention  can  therefore  be  focused  on 
Feshbach  rather  than  orbiting,  resonances.  Theoretical  advances 
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along  these  lines  have  revealed  a  network  of  related  effects,  in¬ 
cluding  transition-state  spectroscopy,  curve  switching  and  line 
narrowing.  Points  of  contact  have  also  been  established  between 
spectroscopy  and  laser-induced  scattering [4] .  In  what  follows,  we 
shall  consider  a  treatment  of  laser-induced  resonances  aimed 
specifically  at  the  problem  of  curve  switching  (in  the  Landau-Zener 
sense).  The  general  problem  will  consist  of  using  a  laser,  LI,  to 
effect  a  transition  from  a  particular  scattering  state  to  a  manifold 
of  bound  states.  A  second  laser,  L2,  is  then  used  to  induce  multi¬ 
photon  transitions  within  this  manifold.  [In  another  context,  L2 
may  be  regarded  as  a  spectroscopic  probe.]  Finally,  the  quasibound 
intermediate  state  decays  by  an  emission  stimulated  by  LI  into  a 
different  channel  (Fig.  1) .  The  reason  for  including  L2  is  that  it 
introduces  an  extra  flexibility  into  the  problem  and  could  lead  to 
considerably  enhanced  probabilities  for  curve  switching,  depending 
on  the  structure  of  the  potential  curves.  This  illustrates  an 
interesting  paradox:  in  the  realm  of  free-to-free  transitions,  the 
introduction  of  an  additional  laser  reflects  a  higher-order  process 
associated  with  a  diminution  of  the  cross-section.  In  the  context 
of  resonance  scattering,  however,  another  laser  can  actually  enhance 
the  cross-section. 

Let  us  consider  the  process  illustrated  in  Fig.  1,  involving 
three  bound  states  |s>,  |t>  and  |u>,  and  two  continua,  |a>  and  |B>, 
dressed  by  two  lasers,  LI  and  L2.  In  the  Breit-Wigner  formulation, 
the  T-matrix  for  curve  switching  is  written  as 


f 


(6) 
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where  Tp  is  the  potential  scattering  term,  arising  from  a 
nonradiative  interaction  (e.g.,  spin-orbit  coupling),  and  TR  is 
the  rapidly  varying  resonant  term  arising  from  the  radiative 
interactions  with  the  manifold  of  bound  states.  It  can  be  shown 
that  T^  can  be  written  as 

Ts«  ■  %  IHPRnHRp'^  >  <7> 

where 


Hpr  *  PHR 

R  a  I  I v><v I 

v=s,t,u 

P  *  l  f  d£*  |  yE'xyE'  | 

a,p  is  a  set  of  angular  momentum  quantum  numbers  describing  the 

open  channels  in  Pig.  1,  s,t,u  is  a  set  of  vibrational  and  angular 

momentum  quantum  numbers  describing  the  three  bound  states  and 

Q  «  lim  [E-H^-H^  (E-H’+iejRpnr1,  (8) 

e+0 


where  H'  -  Hpp  +  HpQ[  (1-R)/ (E-HQQ)  ]HQp,  Q  =  l-p. 

The  matrix  Q"1  can  be  expanded  in  the  basis  spanning  R,  in 
this  case  the  three  bound  states  |s>,|t>  and  |u>.  The  result  is 


where 


rs  -  2*1  | <s |HRR| y> 1 2 

Y“dfp 

Vst  a  <s I Hrr I t> 


(10a) 
(10b) 
etc.  (10c) 


Here#  Eg  is  the  energy  of  the  dressed  state  s. 

It  remains  to  determine  the  in  and  out  scattering  states 

|(l»“>  and  | in  Eq.  (7).  Since  the  matrix  elements  in  (10a)  and 
p  « 

(10c)  are  calculated  in  terms  of  the  standing  wave  representation 
of  the  continuum  states,  it  is  economical  to  use  a  linear 
combination  of  the  same  states  to  describe  l<t>~>  and  |<|,+  >: 

p  « 


1 4>p>  *  I  |y>  e‘i6Y  (11) 

where  is  the  phase  shift  associated  with  channel  y.  The  matrix 
W  can  then  be  determined  from  the  S-matrix  for  potential  scattering 
by 


,2i5y 


W 


y<x 


(12) 


i.e.,  W  is  the  eigenvector  matrix  of  Sp.  Eq.  (7)  can  now  be 
written  in  terms  of  the  W  matrix  as 


v  v 


I 

YY 


Wa  H  Q  • 
'  Py  yv  vv 


H  ,  , 
v  y 


W  ,  e 
Y  « 


i(5y+6y' ) 


(13) 


In  the  distorted-wave  Born  approximation  (DWBA)  there  is  a 
remarkable  simplification  to  Eq.  (13)  since  W  is  the  unit  matrix, 
and  T  may  therefore  be  written  as 
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To  *  £Ho  +  I  Ha  Q  ,  H  ,  “i  e1(6«+6p)  (14) 

pa  S  pa  $v»  pv  vv'  v'a\ 

where  Hpa  is  the  "direct"  or  potential  scattering  terra.  Notice 
that  this  involves  a  definition  of  T  which  is  S  3  l-2xiT.  Another 
definition  common  in  molecular  scattering  is  S  *  1-iT,  in  which  the 
right  hand  side  of  (14)  would  be  multiplied  by  2it.  The  factor  of 

i(6*+6B) 

e  p  is  of  course  irrelevant  to  the  calculation  of  cross- 

sections,  which  depend  upon  |T|2.  it  may  therefore  be  dropped  for 
convenience.  Notice  that  the  EWBA  applies  only  to  potential 
scattering  in  this  case,  and  aside  from  this  TR  is  determined 
exactly. 

To  illustrate  the  possible  scattering  processes  that  occur  in 
Fig.  1.,  we  single  out  three  pathways  giving  rise  to  curve  switching 
in  Fig.  2.  The  total  cross-section  for  switching  from  a  to  p  is 
proportional  to  |T]J2  +  | T2 1 2  +  IT3I2,  where  Ti  is  the  T  matrix  for 
process  i.  The  reason  that  these  three  pathways  are  independent  is 
that  they  give  rise  to  different  final  state  energies.  In  process  1, 
energy  i  is  "conserved"  so  that  there  is  a  potential  scattering 
contribution  which  is  missing  from  processes  2  and  3.  The  term 
Tp  could  arise  from  spin-orbit  coupling  and  is  denoted  by  the 
curved  arrows  in  Fig.  2. 

The  T-matrices  can  be  determined  in  a  consistent  way  simply  by 
inverting  Q”l.  We  find; 


<Tl>Pa 


(Tp) 


pa 


■fi**  i  (v  i  ru)~ v 


2 

tu 


H  H  „ 
ps  as 


(15a) 
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<VPa’ 


'at  (Au  +  5  ru 


H0sH«t 


(15b) 


where 


<T3>Pa 


'tu(*s  +  I  FS, 

D 


VHau 


(*s+5rs)[(Xt+3rt)  (Su+5ru)-v?u]-  v.t(Mrt)  • 


(15c) 


(15d) 


Notice  that  in  each  of  these  terms  the  effect  of  all  the  possible 
ancillary  absorptions  and  emissions  (indicated  by  the  dotted  arrows 
in  Fig.  2 )  are  implicitly  taken  into  account.  It  is  clear  from 


(15)  that  as  the  energy  defects,  A,  of  the  dressed  states  approach 
zero,  the  resonance  states  can  have  a  pronounced  effect  on  the 
magnitude  of  the  cross-section.  [The  energy  defect  in  this  context 
is  essentially  the  detuning  of  L2  from  exact  resonance.]  What  is 
particularly  interesting  is  that  T3  has  the  same  order  of  magnitude 
as  Ti  even  though  it  is  a  four-photon  process!  It  is  possible  that 
under  the  appropriate  circumstances,  T3  might  greatly  exceed  Ti  and 
T2*  Thus,  the  introduction  of  L2  introduces  greater  flexibility  into 
RAIC  and  can  be  used  to  drive  processes  which  are  not  possible  with 
a  single  laser. 

It  is  worthwhile  mentioning  the  advantages  of  a  formal  approach 
to  resonance  scattering,  as  opposed  to  a  coupled-channels  procedure. 


First,  the  radial  portions  of  the  wave  functions  used  to  calculate  the 
matrix  elements  of  (15)  are  independent  of  the  total  angular 
momentum  J.  This  means  that  the  angular  momentum  couplings  can  be 
performed  algebraically.  In  the  coupled-channels  procedure,  the 
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entire  scattering  wave  function  is  propagated  at  once,  so  that  blocks 
of  different  J  are  coupled  together.  For  process  3  in  Fig.  2,  for 
example,  nine  J  blocks  would  be  coupled,  and  obviously  this  would  make 
a  coupled-channels  calculation  both  complicated  and  time  consuming. 
Second,  for  a  realistic  scattering  experiment  there  will  always  be 
a  distribution  of  collision  energies  which  must  be  averaged  over  to 
give  the  observed  cross-section.  To  do  this  numerically,  with  a 
separate  calculation  at  each  collision  energy,  would  clearly  render  a 
coupled-channels  calculation  practically  impossible. 

Finally,  we  note  that  the  effects  of  line  narrowing  [5] , 
previously  discussed  in  reference  to  photo-ionization,  are  also 
evident  in  the  context  of  laser- induced  resonance  scattering.  This 
effect  arises  from  the  interaction  between  (or  "overlapping"  of) 
closely-spaced  bound  states  and  results  in  a  narrowing  of  the  peak  in 
the  scattering  spectrum,  instead  of  the  broadening  which  would 
otherwise  be  expected  because  of  the  interaction  between  the  states. 

III.  Surface  Processes:  Transitions  to  Surface  States  of  Semiconductors 

For  a  finite  semiconductor,  the  energy  levels  are  divided  into 
valence,  conduction  and  surface  bands  [6].  Electrons  in  valence  or 
conduction  states  are  delocalized  throughout  the  solid.  The  surface 
states,  however,  are  confined  to  the  vicinity  of  the  surface. 
Consequently,  this  suggests  the  use  of  a  laser  to  excite  electrons 
from  the  valence  band  into  surface  states  in  order  to  increase  the 
charge  in  the  surface  region.  This  charge  can  provide  a  strong 
Coulombic  interaction  with  an  adspecies. 
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Within  a  simple  one-dimensional  model  and  the  nearly-free- 
electron  approximation,  we  [7]  have  calculated  the  energy 
levels  of  silicon  which  are  given  in  Fig.  3*  The  valence  band  (V) 
is  occupied  in  the  ground  state  whereas  the  conduction  band  (C)  and 
the  surface  band  (S)  are  empty.  The  wave  functions  associated  with 
the  valence  states  are  represented  by  a  sum  of  plane  waves  with 
crystal  momentum  k,  whereas  the  surface-state  wave  functions  are 
damped  oscillators  with  crystal  momentum  g/2  +  i< .  The  probability 
of  a  transition  from  the  valence  state  to  a  surface  state  is  given 
by 

M  »  I <k I  A  *  P  I k> I 2,  (16) 


where  A  is  the  vector  potential  of  the  laser  radiation  and  P  is 
the  momentum  operator  of  the  electron.  Using  the  wave  functions 
in  this  expression,  we  note  that  transitions  that  conserve  the 
real  crystal  momentum  are  favored.  Consequently,  we  can  write 


(17) 


where  I  and  o>  are  the  intensity  and  angular  frequency  of  the  laser 
radiation,  the  subscript  "0"  indicates  integration  over  the  first 
unit  cell  and  z  is  the  coordinate  along  the  crystal  axis.  To 
first  order,  the  absorption  cross-section,  o,  can  now  be  written  as 


/-V/g 

•  *  (!)  (V)f*  Jo  *  "2 


(18) 
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Here,  V  is  the  matrix  element  of  the  crystal  potential  and  g 
the  reciprocal  lattice  vector.  We  have  evaluated  this  expression  for 
silicon  and  the  results  are  depicted  in  Fig.  4.  Near  the  mode  in 
the  figure,  Ik  I  is  greatest  and,  consequently,  the  surface 
density  is  the  greatest.  Since  the  cros Bisection  near  this  node 
is  quite  appreciable,  the  possibility  of  using  a  laser  to  shift 
charge  from  the  bulk  crystal  into  the  surface  region  is 
favorable. 

If  we  now  look  back  at  our  system  before  exposure  to  the 
laser,  the  charge  density  n^(z),  would  be  given  by  [8] 


n0(z) 


-  E(kj  [^(z))2  » 


(19) 


where  4,^ (z)  is  the  valence  wave  function  and  Ep  is  the  Fermi 
energy  with  wave  vector  kp.  However,  with  the  laser  exciting  a 
surface  state  the  charge  density,  n(z),  would  become 

n(z)  »  nQ(z)  +  UkIzH2  (20) 

where  4k (z)  is  the  surface  wave  function.  Since  the  bulk  state 
that  is  being  excited  is  delocalized  throughout  the  semiconductor, 
to  be  exact  we  would  have  to  remove  an  infintesimal  amount  of 
charge  on  the  right  side  of  the  equation.  Eqs.  (19)  and  (20)  have 
been  solved, and  the  results  for  a  typical  surface  state 
[k  ■  —.5| 2V/g | ]  are  presented  in  Fig.  5.  The  dramatic  increase 
in  surface  charge  (dashed  line)  over  the  ground  state  charge 
(solid  line)  produced  by  the  laser  excitation  is  clearly 
indicated. 
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If  a  charge  adspecies  were  in  the  vicinity  of  this  surface, 
the  classical  change  in  the  interaction  potential  produced  by 
this  excited  surface  state  would  be 


5i}  “  ^  J dz  k(z)l2  e 


2  — X  I  2-  2, 


where  Z  is  the  valence  of  the  adspecies  at  position  z i  and  x 
is  the  Thomas-Fermi  screening  parameter  [9}: 


electron 

where  n  is  the  average Adensity  of  the  solid.  Results  of  Eq.  (21) 
are  depicted  for  various  excited  surface  states  in  Fig.  6. 

Even  for  the  surface  states  away  from  the  gap  center 
[ k a . 1 1 2 V/g | ] »  we  can  see  that  the  change  in  potential  is  quite 
large.  Furthermore,  the  decrease  in  charge  above  the  surface 
would  suggest  that  the  screening  parameter,  Eq.  (22),  would  be 
smaller  than  the  bulk  value.  Consequently,  the  interaction  and 
its  range  given  here  should  be  considered  a  lower  limit  of  the 
true  potential. 

We  have  shown  that  not  only  is  it  possible  to  populate 
surface  states  (Fig.  4),  but  the  induced  surface  charge  is  sub¬ 
stantial  (Fig.  5)  and  the  resultant  Coulombic  interaction  can 
have  a  large  effect  on  charged  adspecies  (Fig.  6).  Consequently, 
a  laser  could  be  used  to  enhance  adsorption  or  desorption  by  con¬ 
trolling  the  surface  charge. 
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Figure  6. 


a)  Case  of  two  open  channels  interacting  radiatively 
with  a  manifold  of  bound  states  by  means  of  a  laser, 
LI.  Transitions  within  the  manifold  are  made  by  a 
second  laser,  L2. 

b)  The  same  problem  in  the  "dressed"  representation. 
Here,  |n^>  is  the  number  state  of  the  Fock 
representation  of  laser  i. 

The  various  pathways  leading  to  the  curve-switching 
transition  p«-a.  The  ancillary  absorptions  and 
emissions  of  photons  from  L2  are  shown  as  dotted 
arrows.  These  indirect  pathways  give  rise  to 
interferences  in  the  overall  T-matrix. 

Dispersion  relationship  in  complex  crystal  momentum 
space  (k+iic)  for  a  finite  linear  chain.  The  valence, 
surface  and  conduction  bands  are  labeled  V,  S  and  C, 
respectively. 

Absorption  cross-section  for  surface  states,  o,  in 
A2  versus  the  frequency  of  the  exciting  laser 
radiation. 

Electron  density  distribution  at  the  surface.  The 
solid  line  represents  the  ground  electronic  state, 
and  the  dashed  line  represents  the  system  with  the 
excited  surface  state  K*-.5l2V/g|  in  the  lower 
branch. 

The  magnitude  of  the  surface  interaction  potential 
(in  millihartrees)  at  various  distances  from  the 
surface.  The  solid  line  represents  the  system  with 
excited  state  k—  l2V/g|  ;  the  dashed  line,  ic«-.5l2V/gl 
and  the  dotted  line,  ic— .1  l2V/g| ,  all  in  the  lower 
energy  branch. 
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